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We present a theory to compute the force due to light upon a particle on a dielectric plane by the 
Coupled Dipole Method (CDM). We show that, with this procedure, two equivalent ways of analysis 
are possible, both based on Maxwell's stress tensor. The interest in using this method is that the 
nature and size or shape of the object, can be arbitrary. Even more, the presence of a substrate can 
be incorporated. To validate our theory, we present an analytical expression of the force due to the 
light acting on a particle either in presence, or not, of a surface. The plane wave illuminating the 
sphere can be either propagating or evanescent. Both two and three dimensional calculations are 
studied. 

PACS numbers: 03.50.De, 78.70.-g, 42.50.Vk, 41.20.-q 
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I. INTRODUCTION 



The demonstration of mechanically acting upon small 
particles with radiation pressure was done by Ashkin and 
coworkers 1,2 . A consequence of these works was the in- 
vention of the optical tweezer for non destructive manip- 
ulation of suspended particles 3 or molecules and other 
biological objects 4-6 . Recently, these studies have been 
extended to the nanometer scale 7-12 , and multiple par- 
ticle configurations based on optical binding have been 
studied 13-17 . Also, the effect of evanescent waves cre- 
ated by total internal reflection on a dielectric surface 
on which particles are deposited was studied in Ref. [18]. 
However, the only theoretical interpretation of such sys- 
tem is given in Refs. [19-20]. In Ref. [19] no multiple 
interaction of the light between the particles and the di- 
electric surface was taken into account. On the other 
hand, in Ref. [20] a multiple scattering numerical method 
was put forward limited to a 2-D configuration. 

It is worth remarking here that several previous theo- 
retical works on optical forces usually employ approxima- 
tions depending on the radius of the particle; if the parti- 
cle is small it has been usual to split the force into three 
parts: the gradient, scattering, and absorbing forces 21 . 
However, a rigorous and exact calculation requires the 
use of Maxwell's stress tensor. We shall use it in this 
paper. Some work has been done in free space 7 ' 22 , or for 
a spherical particle over a dielectric surface illuminated 
by a Gaussian beam. 23 

We shall present, therefore, a detailed theoretical anal- 
ysis in three dimensions of how the optical force is built 
on the multiple interaction of light with the particle and 
the dielectric surface. This will be done whatever its size, 
shape, or permittivity. To this end, we shall make use of 
the Coupled Dipole Method (CDM), whose validity was 
studied in detail in Ref. [24]. 

In Section II we present the CDM, and two possibili- 
ties that arise with this method to compute the force by 



means of Maxwell's stress tensor. Concerning the first 
one, in Section II A we use directly Maxwell's stress ten- 
sor and perform the surface integrations. As regards the 
second one, we present in Section II B the dipole approxi- 
mation on each subunit of discretization for the numerical 
calculations. Since however these methods are somewhat 
cumbersome from a numerical point of view, we have in- 
troduced in Section III an analytical calculation for the 
force due to the light on a small particle in the presence 
of the surface. Results are illustrated in three dimensions 
in Section III A (a little sphere) and in two dimensions in 
Section IIIB (a small cylinder). In Section IV we com- 
pute the force with the CDM, and we validate these cal- 
culations on electrically small particles by means of the 
analytical solution presented in Section III. After this 
validation of the CDM on little particles, we present in 
Section IV C calculations on larger particles. 



II. ELECTROMAGNETIC FORCE COMPUTED 
WITH THE COUPLED DIPOLE METHOD 

The Coupled Dipole Method (CDM) was introduced 
by Purccll and Pennypacker in 1973 for studying the 
scattering of light by non-spherical dielectric grains in 
free space. 25 This system is represented by a cubic array 
of N polarizable subunits. The electric field E(rj,a;) at 
each subunit position can be expressed as: 



E(r-i, a) = E (r i} w) + V[S(r i; rj , oj) 



N 



(i) 



where Eo (r^ , ui) is the field at the position in the ab- 
sence of the scattering object, T is the field susceptibility 
associated to the free space 26 (with T(r^,r$,a;) = to 
avoid the diagonal case), and S represents the field sus- 
ceptibility associated with the surface in front of which 
the particle is placed (see Fig. 1). The derivation of the 
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field susceptibility of the surface is extensively developed 
in Refs. [ 27, 28]. ctj(u)), the polarizability of the subunit 
j, is expressed as: 



,-M=a°( W )/[l-(2/3)iA3a$M] 



(2) 



where fco = |ko| = lo/c (ko being the incident wavevec 

3 



tor of the electromagnetic field in vacuum) and orAu)) is 



given by the Clausius-Mossotti relation : 
3d 3 e{uj) - 1 



4tt £(w) + 2' 



(3) 



In Eq. (3) d is the spacing of lattice discretization and 
e(u) stands for the relative permittivity of the object. 
Let us remark that the polarizability is expressed ac- 
cording to Eq. (2) as defined by Draine 29 . The term 
(2/3)ifcga°(w) is the radiative reaction term, necessary 
for the optical theorem to be satisfied and for a correct 
calculation of forces via the CDM. 30 



A. Force computed with Maxwell's stress tensor 

The force F on an object due to the electromagnetic 
field 31 is computed from Maxwell's stress tensor: 32 



l/(87r)3?e 



[(E(r,w).n)E*(r,w) 



(5) 



(H(r,w).n)H*(r,w) 
l/2(|E(r, W )| 2 + |H(r,c)| 2 )n]dr 



where S is a surface enclosing the object, n is the local 
outward unit normal, * denotes the complex conjugate, 
and 3?e represents the real part of a complex number. Let 
us notice that Eq. (5) is written in CGS units for an ob- 
ject in vacuum, and so will be given all forces presented 
in Section IV. To apply Eq. (5) with the CDM, we must 
first solve Eq. (1) to obtain E(r,,u;) at each dipole po- 
sition, and then, through Eq. (4) and Faraday equation, 
the electromagnetic field is computed at any position r 
of S. This enables us to numerically perform the two 
dimensional quadrature involved in Eq. (5). 




FIG. 1. Geometry of the configuration considered in this 
paper: sphere, or cylinder, of radius a on a dielectric flat sur- 
face. The relative permittivity is e — 2.25 both for the sphere 
(or the cylinder) and the surface. The wavelength used is 
A = 632. 8nm in vacuum and the incident wave vector k is in 
the XZ plane. 

Once the values of E(r^, u>) are obtained by solving the 
linear system, Eq. (1), (whose size is 3N x 3N), it is easy 
to compute the field at an arbitrary position r: 

N 

E(r,w) = E (r,u;) + ^[S(r,r j ,w) (4) 

i=i 

+ T(r,r J ,w)]a J (w)E(r J ,w). 

The computation of the force also requires the magnetic 
field radiated by the scattering object. We obtain it 
through Faraday's equation, H(r, u) = c/ (ioj) V x E(r, u>) 



B. Force determined via the dipolar approximation 

Let us consider a small spherical particle with a radius 
smaller than the wavelength. Then the w-component of 
the force can be written in the dipole approximation: 33 ' 30 

F u (r a ) = (l/2)3ee]T L(r , w) " } ) , (6) 

(«=1, 2, 3) 

where r is the position of the center of the sphere and 
u and v stand for the components along either x, y or 
z. We discretize the object into N small dipoles p»(r, uj) 
(i=l,...,N) so that it is possible to compute the force on 
each dipole from Eq. (6). Hence, to obtain the total force 
on the particle it suffices to sum the contributions F(r^) 
from of all the dipoles. To use this method it is necessary 

to know d^'"^ 1 ' 4 ' ^ a t each discretization subunit. On 
ou 

performing the derivative of Eq. (1) we obtain: 

N ( d \ 
+ Ek [S(r ' rj ' w) + T(r,r ^ w)1 a»E(r j; a;). 

j=l V /r=r 4 

Thus, the derivative of the field at r iy requires that 
of E (ri,cj) and that of the field susceptibility both in 
free space and in the presence of the surface for all pairs 
(rj,rj). Hence we now have two tensors with 27 compo- 
nents each. It is important to notice that the derivative 
of the field at Vi has been directly computed from just the 
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field at this position rj, so it is not computed in a self- 
consistent manner. To have the required self-consistence 
for the derivative, it is necessary to perform in Eq. (1) 
a multipole expansion up to second order. Then, this 
equation must be written up to the quadrupole order af- 
ter taking its derivative. As a result, we obtain a linear 
system, whose unknowns are both the electric field and 
its derivative. The disadvantage of this method is that 
the size of the linear system increases up to 12 N x 12N 
and requires the computation of the second derivative of 
the field susceptibility (81 components). More informa- 
tion about the CDM by using the multipole expansion 
can be found in Ref. [24]. 

In what follows, we shall denote CDM-A the force com- 
puted directly from Maxwell's stress tensor Eq. (5) and 
CDM-B the force obtained on using the field derivative 
Eq. (6). The advantages of these two methods is that 
they are not restricted to a particular shape of the ob- 
ject to be discretized. Furthermore, this object can be 
inhomogcncous, metallic, or in a complex system when- 
ever it is possible to compute its field susceptibility. 



. 2 3Ffe / 4zQa(ik z 



"I- EC) y I I 

1 0zl 2 V 4sg 
for p-polarization and 



ctA 



6zg|a| 2 A 
|4z 3 + aA| 2 



F x 



3?e 



2 



lz$a(ik x )* 
xA 



8zl 



F z = \E , 



3?e f8z$a(ik z y 



2 V 8z o + aA 



+ 



12zi 



2 A 



8z3 + aA| 2 



(11) 
(12) 



for s-polarization. We see that the advantage of working 
with the static approximation is that an analytic form of 
the force is obtained. To see the effect of the incident 
field only (i.e., without interaction with the surface), we 
can put z — > oo or A = in Eqs. (9)-(12). The forces 
are then expressed as: 



E \ 2 ^f (a(ik x y). 



(13) 
(14) 



III. FORCE ON A DIPOLAR PARTICLE 

A. The three dimensional case: a sphere 

Eq. (1) with N = 1, taking the surface into account, 
gives for the field at the position r = (xo,yo,^o) of the 
sphere of a radius a: 



\Eo y \ 2 for s-polarization and 



\En\ 2 = 



E(r , w) = [I - a(w)S(r , r , w)] 1 E (r , uj), 



(8) 



where I is the unit tensor, and a(uj) the polarizabil- 
ity of the sphere according to Eq. (2) with ao(co) = 
a 3 (e(uj) — 1)/(e(lj) + 2). We notice that S is purely di- 
agonal and depends only on the distance zq between the 
center of the sphere and the surface (see Fig. 1). We also 
assume that the sphere is near the surface and, hence, the 
field susceptibility of the surface can be used in the static 
approximation (fco = 0, we shall discuss the validity of 
this approximation in Section IV). Therefore, the com- 
ponents of this tensor become S xx = S yy = —A/(8>Zq), 
and S zz = -A/(4z$), with A = (1 - e)/(l + e) repre- 
senting the Fresnel coefficient of the surface. Since we 
consider the object in the presence of a surface with a 
real relative permittivity, A is real. As shown by Fig. 1, 
the light incident wave vector k lies in the XZ plane. 
Therefore, there is no force in the Indirection. On us- 
ing Eqs. (6) and (8), and assuming the incident field Eo 
above the surface to be a plane wave either propagating 
or evanescent, depending on the illumination angle 9, the 
components of the force on the sphere can be written as: 



F x = 



9fJe 
~~2 



1az%{ik x 



2\Eo 



\En 



izX + aA 



4zl 



aA 



F, = \E t 



Oxl 



3?e (%zla{ik z )* 12z, 



+ 



2 \8z§ + aA |8zg + aA 



(9) 
(10) 



with \E \ 2 

\E 0x \ 2 + \E 0z \ 2 for p-polarization. Eqs. (13)-(14) show 
a spherical symmetry, and hence the results both in p 
and s-polarization are the same. 

If we look at Fig. 1, we see that the incident field above 
the surface always has k x real, but k z can be either real 
(propagating wave) or imaginary (evanescent wave when 
9 > 9 C where 9 C is the critical angle defined as y/e sin 9 C — 
1). Hence, all forces in the A-direction have the form 
Adie(a(ik x )*) where A is always a positive number. In 
using Eq. (2) we find that ^te(a(ik x )*) ~ (2/3)alk%k x 
(we have assumed that {A/^k^a^ <C 1, in fact this ex- 
pression is about 6.6 x 10~ 7 for a = lOnm, A — 632. 8nm, 
and e = 2.25, thus this approximation is perfectly valid). 
Hence, whatever the field, either propagating or evanes- 
cent, and the system either in presence of a surface or in 
free space, the force in the A-direction is always along 
the incident field. 

From Eq. (14) and from the discussion above, it is easy 
to see that in the absence of interface the force is pos- 
itive for a propagating incident wave (k z real). In the 
case of an evanescent incident wave, k z = ij with 7 > 0, 
and hence the force becomes F z = — jao\Eo\ 2 /2; namely 
the sphere is attracted towards the higher intensity field. 
Concerning the force along the Z-dircction, its sign will 
depend on the nature of the field and the interaction of 
the sphere with the surface. We shall discuss this in Sec- 
tion IV A. 



B. The two dimensional case: a cylinder 

For a cylinder with its axis at (xq, zq), parallel to the 
F-axis (Fig. 1), the electric field at its center is obtained 
by an equation similar to Eq. (8), but with a different 
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polarizability. With the help of Refs. [34] and [35] we 
write this polarizability: 



ai(w) 

Ct2(uj) 



ifcg7ra;(w)/2' 
a§(w) 



with a?(w) 



e(a;) - 1 a 2 

+ 1 2 



(15) 



as those established for the sphere with only a replace- 
ment of a by oil or a 2i depending on the polarization. 
Concerning the force along the AT-dircction, we have the 
same effect as for the sphere, namely F x has the sign of 
k x • 



l-ifcg7ra§(w)' 



with a° 2 (Lu) = (e(w) - 1) — . (16) 



The subscripts i = 1, and 2 correspond to the field ei- 
ther perpendicular or parallel to the axis of the cylin- 
der, respectively. The field susceptibility of the sur- 
face in the two dimensional case is given in Ref. [36] 
for s-polarization and in Ref. [37] for p-polarization. 
Since we address a cylinder with a small radius a and 
near the surface, we use the static approximation, then 
S X x = S zz = — A/(2zg) and S yy = 0. In the same way as 
seen before, the force is written as: 



F x = \E \ 2 ^(a 2 (ik x y) 
F t = \Eo\ 2 ™ (a 2 (ik z y). 



for s-polarization and 



Fx = \Er\ 



F 7 



»e (2z$a 1 (ik x )* 
2z 2 + aiA 
2z 2 a 1 (ik z y 



2 



2^o|ai| 2 A 



2 V 2z£ + ai A |2z^+aiA| 



(17) 
(18) 

(19) 
(20) 



for p-polarization. \E$\ — \Eq v \ for s-polarization and 
|_E | 2 = \Ev x \ 2 + l-^Ozl 2 for p-polarization. If, again, 
z — > oo or A = 0, there is no interaction between the 
cylinder and the surface, then we find the same equations 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this section we present numerical results on forces 
acting on either a small sphere or a small cylinder. The- 
ses forces are normalized in the form F U /\E \ 2 where F u 
is the w-component of the force, and \E \ stands for the 
modulus of the incident field at the center of either the 
sphere or the cylinder. All calculations are done for a 
body in glass (e — 2.25), at a wavelength of 632. 8nm, in 
front of a flat surface (e = e = 2.25) illuminated from the 
glass side by internal reflection (Fig. 1). 



A. Results for a small sphere 

We have first checked our CDM calculation by com- 
paring it with the well known Mie scattering results for a 
sphere in free space illuminated by a plane wave. 38 The 
force is: 



FMie — —\Eo\ 

87T 



COS OCsca) ~T 



(21) 



where C ex t denotes the extinction cross section, C sca the 
scattering cross section, and cos the average of the co- 
sine of the scattering angle. Calculations are done for a 
sphere of radius a — lOnm. 



CDM-A 


; CDM-B 


j dip. app. 


Mie 


force 


N 


%(Mie) 


force 


N 


%(Mie) 


force 


%(Mie) 


force 


2.8119E-22 


81 


0.46 


2.8338E-22 


81 


1.24 


2.8027E-22 


0.13 


2.7991E-22 


2.8181E-22 


912 


0.68 


2.8243E-22 


912 


0.91 


2.8151E-22 


1791 


0.57 


2.8194E-22 


1791 


0.73 


2.8151E-22 


2553 


0.57 


2.8186E-22 


2553 


0.70 



TABLE I. Force on a sphere of radius a = lOnm in free space. Numerical results for different number of subunits N in 
the CDM-A, CDM-B. Comparison of calculation with the dipolar approximation and Mie's calculation. %(Mie) is the relative 
difference (in percent) between the exact Mie calculation and the method used. 



Table I compares the force obtained from the CDM 
on using, without any approximation, either the method 
developed in Section II A (CDM-A) or that from Sec- 
tion II B (CDM-B), and from the dipolar approxima- 
tion (dip. app.) presented in Section III A, with the 
Mie calculation (%(Mie) is the relative difference in per- 
cent between the Mie result and the other corresponding 
method). For an incident field with \E \ = 94825V/m, 



which corresponds to a power of 1.19mW distributed on 
a surface of lO^im 2 , the force on the sphere in MKSA 
units is 2.7991 x 10~ 22 Newtons. One can see that for 
both CDM-A and B the convergence is reached even for 
a coarse discretization, and, hence, either one of the two 
CDM approaches can be used. As regards the dipolar 
approximation, we conclude that it is perfectly valid to 
use it for a sphere of radius a =10nm (a/A < 0.016). 
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Now that we have validated our methods (both ana- 
lytic and CDM) we proceed to take the surface into ac- 
count. It should be remarked that with the CMD-A it 
is not possible to compute the force when the sphere is 
on the surface. This is because for an observation point 
very close to the sphere, the electromagnetic field values 
are affected by the discretization of the sphere, and so 
the field is not correctly computed. An empirical crite- 
rion that we have found 24 is that the electric field must 
be computed at least at a distance d from the sphere 
but this criterion depends on the relative permittivity. 
For more precision about the dependence of the criterion 
and the relative permittivity one can look to Ref. [29]. 
With the CDM-B this problem does not occur because 
with this approach it is not necessary to obtain the field 
outside the sphere. 



a sphere of radius a = lOnm, at an angle of incidence 
9 = 42°, without any approximation with the CDM- A 
(namely, taking into account all retardation effects) with 
the SAFSS, and with the approximation in which no in- 
teraction between the sphere and the surface is consid- 
ered. The difference between SAFSS and the exact cal- 
culation is less than 1.5%. This is in fact logical. Near 
the surface, the SAFSS is correct, far from the surface, 
however, the field susceptibility associated with the sur- 
face in the exact calculation is significantly different from 
the field susceptibility derived from a static approxima- 
tion. Nevertheless, for distances larger than z = 30nm 
the curves overlap because the sphere does not feel the 
substrate at this distance. This is manifested by a differ- 
ence of only 2% between the exact calculation result and 
that computed without addressing the surface (horizon- 
tal line). 




-10 



10 20 30 40 50 
distance between the sphere and the surface (nm) 



FIG. 2. Normalized force in the Z-direction on the sphere 
of a — lOnm versus distance Z. The angle of incidence of 
illumination is 6 = 42° in p-polarization. The full line rep- 
resents the exact calculation with CDM-B, the dashed line 
corresponds to the static approximation with CDM-B, and 
the dotted line is the calculation without interaction between 
the sphere and the surface. 

In all figures shown next, we plot the force versus the 
distance z between the sphere (or the cylinder, see Sec- 
tion IV B) and the plane, (notice that we represent by z 
the distance between the center of the sphere, or cylin- 
der, and the plane). The calculation using the dipole 
approximation, as well as the CDM(A-B), has been done 
with the static approximation for the field susceptibility 
of the surface (SAFSS). However, the distance between 
the sphere and the surface goes generally up to lOOnm. 
In order to justify the study of the force at distances 
about lOOnm between the sphere and the plane through 
a calculation done in the static approximation, we plot 
in Fig. 2 the normalized force F z for p-polarization, with 




Mie limit 



zero force /47nm 



50 100 

distance between the sphere and the surface (nri) 



-5 

25 50 75 100 

distance between the sphere and the surface (nm) 



FIG. 3. Normalized force in the ^-direction on a sphere of 
radius a = lOnm. The full line corresponds to the dipole ap- 
proximation, the dashed line to the CDM-A, and the dotted 
line to the CDM-B. The angle of incidence is 6 = 0°. The 
inset shows the force near z = 50nm. We show the zero force 
and the force computed from Mie's limit with Eq. (21). 

Fig. 3 shows the normalized force for light at an an- 
gle of incidence 9 = 0°. The curves corresponding to 
CDM-A and B are similar, and the dipole approximation 
appears slightly above when the sphere is close to the 
surface. This may seem strange at first sight in view of 
the good results presented in Table I (we will discuss it 
later). We can see that although the illuminating wave 
is propagating, if the sphere is near the surface, it is 
attracted towards it, opposite to the propagation direc- 
tion. To understand this, we look at Eq. (12), established 
with the dipole approximation with the values k x = 0, 
k z = k , E 0z = which corresponds to 9 = 0°. Af- 
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ter some approximations (namely (4/9)^0^ <C 1 which 
implies |a| 2 ~ a 2 ,), the force can be written: 



|£p| 2 64zg 
1 8^ 



aA\ 



32z,f 



(22) 



The factor before the bracket of Eq. (22) corresponds to 
the intensity of the field at the position of the sphere. 
The first term in the bracket of this equation is due to 
the light scattering on the particle (as in free space) and 
is always positive. The second term in the bracket is al- 
ways negative as A < 0. Therefore, the relative weight 
of the two terms in Eq. (22) determines the direction of 
F z . F z given by Eq. (22), becomes zero for: 



9(g~l) 
32^( £ + l)- 



(23) 



Hence, in our example we find z = 57nm. Below the 
value of Eq. (23) the force is attractive towards the sur- 
face, and above this value the sphere is pushed away. 
This is seen in the inset of Fig. 3 which enlarges those 
details. We find F z — at z = 47nm namely at 
zq =(47+10)nm=57nm, which is exactly the same value 
previously found. Physically, the attraction of the sphere 
is due to the second term of Eq. (22) which corresponds 
to the interaction of the dipole with its own evanescent 
field reflected by the surface. Now we can explain the 
discrepancy between the dipole approximation and the 
CDM as regards the good results obtained in free space. 
In fact, when the computation is done in free space the 
field can be considered uniform over a range of 20nm. 
However, in an evanescent field, the applied field is not 
uniform inside the sphere and the Clausius Mossotti re- 
lation is less adequate, hence, the dipole approximation 
departs more from the exact calculation. However, when 
the sphere is out from the near field zone, the three meth- 
ods match together (see the inset of Fig. 3). We can also 
see in the inset of Fig. 3, that these three curves tend 
towards the Mie limit because at large distance there is 
no interaction with the surface. 




-10 



10 20 30 

distance between the surface and the sphere (nm) 



FIG. 4. Normalized force in the Z-direction acting on the 
sphere with a =10nm. The angle of incidence 6 = 42° is 
larger than the critical angle 9 C = 41.8°. The full line cor- 
responds to the dipole approximation, the dashed line to the 
CDM-A, and the dotted line to the CDM-B. Curves without 
symbols are for p-polarization, and those with symbol + arc 
for s-polarization. 

Fig. 4 shows the ^-component of the normalized force 
when the incident wave illuminated at 9 = 42° > 41.8° = 
6 C . Then, for s-polarization we can write Eq. (12) as: 



F z 



\En 



|8z3 + aA| 2 



[-4*g7ao(a A + 8z 3 ) + 6«gagA] . (24) 



It is easy to see that for a dielectric sphere both the first 
and second terms within the brackets of Eq. (24) are 
always negative. Hence, the sphere is always attracted 
towards the surface (the same reasoning can be done for 
p-polarization). Near the surface the force becomes larger 
because of the interaction of the sphere with its own 
evanescent field. We notice that the normalized force 
becomes constant at larger z. This constant reflects the 
fact that the force decreases as e~ 27Z from the surface. 



B. Results for a small cylinder 

Let us now address an infinite cylinder. Since the CDM 
method used here works in three dimensions, we have 
computed the force on a finite length cylinder. In order 
to verify this approximation, we once again compare the 
force, obtained in free space from the CDM with different 
cylinders lengths, with that from a calculation done with 
the dipole approximation established in Section IIIB, 
and that from an exact calculation for an infinite cylin- 
der 38 (i.e. the well known 2-D version for cylinders of 
the Mie calculation for spheres, hence referred to in the 
table as "Mie") . We consider a radius of the cylinder, 
a = lOnm, with the same spacing lattice as for the case of 
the sphere, namely 81 subunits. We have seen that this 
value of d = 4nm gives consistent results. In all cases we 
compute the force per unit length of the cylinder. 

The first case addressed is with the electric field per- 
pendicular to the axis of the cylinder (p-polarization) . 
The results are given in Table II. The second case con- 
sidered (s polarization) in Table III. 

We notice that the dipole approximation gives cor- 
rect results for p-polarization, but it is worse for s- 
polarization. If we compare the CDM-A and CDM-B, 
we see that they both give the same results. But we 
also see that the length of the cylinder has a great in- 
fluence on them, although up to a different extent ac- 
cording to whether we deal with p or s-polarization. For 
p-polarization, the simulation of an infinite cylinder be- 
comes correct at L ~ A/2 and for s-polarization only it 
is so at L ~ 2A. This can be understood by the fact 
that in p-polarization the electric field is continuous at 
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the end of the cylinder, thus, the end has not a large in- 
fluence on the field computed around the cylinder. How- 
ever, in s-polarization the field is discontinuous at the 
end of the cylinder and then the field will strongly vary 
around this end and so will do the force. This is why in 
s-polarization it is necessary to consider cylinders with 



large lengths in order to avoid edge effects. Now let us 
address the presence of the plane surface to compute the 
force. We consider the cylinder length L = 1551nm. Like 
for the sphere, we address both 6 = 0°, (Fig. 5), and 42°, 
(Fig. 6). The curves from CDM-B stop at z =10nm due 
to the disadvantage previously remarked. 



CDM-A 


| CDM-B 


! dip. app. 


Mie 


force 


L(nm) 


%(Mie) 


force 


L(nm) 


%(Mie) 


force 


%(Mie) 


force 


2.1540E-13 


197 


24 


2.1625E-13 


197 


24 


2.8433E-13 


0.27 


2.8354E-13 


2.9906E-13 


391 


5.47 


3.0013E-13 


391 


5.85 


2.7907E-13 


777 


1.58 


2.8000E-13 


777 


1.25 


2.8661E-13 


1164 


1.08 


2.8756E-13 


1164 


1.42 


2.8347E-13 


1551 


0.03 


2.8439E-13 


1551 


0.30 



TABLE II. Force on a finite cylinder of radius a = lOnm in free space. The discretization interval is d = 4nm. Numerical 
results are presented for different lengths L of the cylinder for both CDM-A and CDM-B. Comparison is made with both the 
dipolar approximation and Mie's calculation. %(Mie) is the relative difference between the exact Mie calculation for an infinite 
cylinder and the method used. Calculations are done for the field perpendicular to the axis of the cylinder. 



The second case considered is with the electric field parallel to the axis of the cylinder (s-polarization): 



CDM-A 


| CDM-B 


; dip. app. 


Mie 


force 


L(nm) 


%(Mie) 


force 


L(nm) 


%(Mie) 


force 


%(Mie) 


force 


0.5649E-12 


197 


63 


0.2163E-12 


197 


86 


1.5015E-12 


2.31 


1.5370E-12 


0.9986E-12 


391 


35 


1.0021E-12 


391 


35 


1.3059E-12 


777 


15.0 


1.3103E-12 


777 


14.7 


1.3971E-12 


1164 


9.10 


1.4018E-12 


1164 


8.80 


1.4430E-12 


1551 


6.12 


1.4479E-12 


1551 


5.80 



TABLE III. The same as in Table II but for the electric field parallel to the axis of the cylinder. 



o 



s polarization 



■p polarization 



-5 

25 50 75 100 

distance between the cylinder and the surface (nm) 



FIG. 5. Normalized force in the Z-direction on a cylinder 
with radius a =10nm, A =632. 8nm, and e = 2.25. The light 
angle of incidence is 6 = 0°. The full line corresponds to the 
dipole approximation, the dashed line to the CDM-A, and the 
dotted line to the CDM-B. 



Concerning Fig. 5, if we focus on F z for p-polarization, 
we can write this force approximated from Eq. (20) by: 



F, = 



2z% 



QlAI 



(a°) 2 fc V4 + 



K) 2 A 

4z 3 



(25) 



Eq. (25) is of the same form as Eq. (22). Hence, the same 
consequence is derived: near the surface the cylinder is 
attracted towards the plane surface. But far from the 
plane the cylinder is pushed away because at this dis- 
tance the cylinder cannot interact with itself. Like for 
the sphere, we can compute the distance zq at which the 
force is null: 



(e-l) 



(26) 



which in our illustration leads to z = 50nm. Although 
we do not present now an enlargement with details of 
Fig. 5, we have found the value zo=(40+10)nm=50nm. 
The cylinder in p-polarization has the same behaviour as 
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the sphere. However, in s-polarization there is a differ- 
ence. Then the force obtained from the dipolar approxi- 
mation is always constant because there is no interaction 
with the surface. This is clear from Eq. (19), and it is due 
to the fact that in the electrostatic limit the field suscep- 
tibility S yy tends to zero, and then there is no influence of 
the surface on the cylinder. This is a consequence of the 
continuity of both the field and its derivative of both the 
plane and the cylinder. 36 Therefore, the cylinder does 
not feel the presence of the plane, and as the wave is 
propagating, the force is positive thus pushing the cylin- 
der away from the plane with magnitude values given by 
table III. Notice that the force obtained from CDM-B, 
when the sphere is in contact with the surface, becomes 
negative in s-polarization. This is due to the diffraction 
of the field at the end of the cylinder, which induces a 
component perpendicular to the plane, and therefore an 
attractive force. 



p polarization 




s polarization 



25 50 75 100 

distance between the cylinder and the surface (nm) 

FIG. 6. Normalized force in the Z-direction on the same 
cylinder as describes in Fig. 5 but with an angle of incidence 
9 = 42° larger than the critical angle C = 41.8°. The full line 
corresponds to the dipole approximation, the dashed line to 
the CDM-A, and the dotted line to the CDM-B. 

In the case represented in Fig. 6, like for the sphere, we 



observe a force always attractive (F z < 0) whatever the 
polarization. For p-polarization we have exactly the same 
behaviour as for the sphere. However, for s-polarization 
the normalized force is always constant whatever the dis- 
tance between the cylinder and the surface, due to the 
same reason as before, namely, S yy = 0. Only when the 
cylinder is on the surface, we can see from the CDM-B 
calculation that the force is slightly more attractive for 
the same reason previously quoted. 



C. Results for a sphere beyond the Rayleigh regime 

Let us now consider a sphere of radius a =100nm. This 
size is far from the Rayleigh scattering regime (~ A/3). 
As in previous cases, we first validate our method with 
the aid of Mie's calculation in free space. The following 
table shows the results: 

As before, as d decreases, the CDM results tend to the 
Mie's calculation. The error never exceeding 1.7%. Now, 
we address the presence of a flat dielectric surface. The 
forces, to be shown next, are computed with CDM-B only 
since the particle can be in contact with the surface. 
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FIG. 7. Normalized force in the ^-direction on a sphere 
with radius a =100nm, A =632. 8nm, and e = 2.25. The light 
angle of incidence is 6 = 0°. The full line corresponds to the 
exact calculation with CDM-B, and the dashed line represents 
the static approximation. 



CDM-A 


CDM-B 


Mie 


force 


N (d =nm) 


%(Mie) 


force 


N (d =nm) 


%(Mie) 


force 


2.1355E-16 


280 (25) 


1.31 


2.1439E-16 


280 (25) 


1.71 


2.1080E-16 


2.1353E-16 


912 (17) 


1.30 


2.1402E-16 


912 (17) 


1.53 


2.1332E-16 


1791 (13) 


1.20 


2.1367E-16 


1791 (13) 


1.37 


2.1312E-16 


4224 (10) 


1.11 


2.1333E-16 


4224 (10) 


1.21 



TABLE IV. Force on a sphere of radius a — lOOnm in free space. Numerical results are for different number of subunits N 
in the CDM-A, CDM-B. Comparison with Mie's calculation. 
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In Fig. 7 we present the case 8 — 0°. We have plotted 
two curves: the exact calculation and the SAFSS done 
with N = 1791. In the inset of Fig. 7, we see that even 
near the surface the SAFSS is not good. This is due 
to the large radius of the sphere, then the discretization 
subunits on the top of the sphere are at lOOnm from the 
surface, and thus the effects of retardation are now im- 
portant. The SAFSS calculation also shows that at a 
distance of 200nm (which corresponds to the size of the 
sphere: 2a=200nm) the sphere does not feel the surface, 
as manifested by the fact than then the curve obtained 
from this computation reaches the Mie scattering limit 
previously obtained in Table IV (cf. the full horizon- 
tal line in the inset). Hence, we conclude that evanes- 
cent waves are absent from the interaction process at 
distances beyond this limit. From the exact calculation 
we obtain a very low force near the surface, due to the 
interaction of the sphere with itself. This effect vanishes 
beyond z w 50nm where oscillations of the force F z take 
place with period A/2. As these oscillations do not occur 
in the SAFSS, this means that they are due to interfer- 
ences from multiple reflections between the surface and 
the sphere. As expected, they decrease as the sphere goes 
far from the surface. 
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FIG. 8. Normalized force in the Z-direction on a sphere 
with radius a =100nm, A =632. 8nm, and e = 2.25. The light 
angle of incidence is 9 — 42° > 6 C . The full line corresponds 
to the exact calculation with CDM-B, and the dashed line to 
the static approximation, the curves without symbol are in 
p-polarization, and those with the + symbol in s-polarization. 

Fig. 8 shows the force computed with an angle of inci- 
dence 8 = 42°. We plot the exact calculation (full line) 
and the SAFSS (dashed line) both for p-polarization (no 
symbol) and s-polarization (+ symbol). Once again, we 
see that the SAFSS is not adequate even near the sur- 
face. On the other hand, in the exact calculation, the 
two polarizations show oscillations of the force F z with 
period A/2. However, there is a large difference of mag- 
nitude of these oscillations between the two polarizations 
(see inset of Fig. 8). To understand this difference, we 
must recall that the sphere is a set of dipoles. When a 
dipole is along the Z-direction there is no propagating 
wave in this direction. But if the dipole is oriented in 



the X(or V)-direction, its radiation is maximum in the 
Z-direction. However, in s-polarization all dipoles are, 
approximately, parallel to the surface, so there is an im- 
portant radiation from the dipole in the Z-direction, and, 
consequently, between the sphere and the surface. 



V. CONCLUSIONS 



In this paper we have presented exact three dimen- 
sional calculations based on the Coupled Dipole Method 
and an analytical expression for the force on either a 
sphere or an infinite cylinder, both in front of a flat di- 
electric surface. The results for small bodies show that, 
whatever the polarization, in the case of a sphere, and 
in p-polarization for the cylinder, the force always has 
the same behaviour: namely, in the case of illumination 
under total internal reflection, the particle is always at- 
tracted towards the surface. A surprising result in the 
case when the illuminating beam is perpendicular to the 
surface and the object remains sticked to the surface, is 
that then the force is attractive due the interaction of the 
particle with itself, and therefore this object keeps sticked 
to the surface. However, when the object is far from the 
surface, the force becomes repulsive, as one would have 
expected. 

For s-polarization, the cylinder does not "feel" the 
presence of the substrate. This is more noticeable for 
a propagating wave, namely, at angles of incidence lower 
than the critical angle. However, when an evanescent 
wave is created by total internal reflection, the force is 
attractive under s-polarization. 

The scope of the static calculation for this configura- 
tion has been validated. We have also shown the ad- 
vantage of having an analytical form which shows the 
contribution of the incident field on the particle, as well 
as that of the force induced by the sphere (or cylinder) on 
itself, thus yielding a better understanding of the physical 
process involved. 

For bigger spheres, we have observed somewhat dif- 
ferent effects of the forces. Under the action of evanes- 
cent waves, the force is always attractive, but it always 
becomes repulsive when it is due to propagating waves. 
Unlike the case of small sphere, there is no point of zero 
force. 
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